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Abstract This article addresses the splitting failure of 
moment-resisting dowel-type fastener joints, in which the 
failure may be attributed to the perpendicular-to-grain 
loading of one single dowel located close to the end of a 
beam. A quasi-non-linear fracture mechanics model based 
on beam on elastic foundation theory is applied. A simple 
approximation suitable for practical design is also proposed. 
Model predictions of the influence of edge distances and 
end distances are compared with test results. 
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Introduction 

Dowel-type fastener joints fail either in a ductile manner, 
characterized by bending of the fastener and/or embedment 
of the fastener into the wood, or in a brittle manner, charac¬ 
terized by splitting of the wood. The ductile failure modes 
are fairly well understood. Widely accepted design expres¬ 
sions have been derived for individual fasteners based on 
theory of plasticity, 1,2 and also the strength and stiffness of 
moment-resisting joints can be reliably predicted. 3,4 Re¬ 
search on the brittle splitting failure modes is still at an early 
stage. The splitting failure mode of centrally loaded joints 
has only recently been studied analytically, 5 9 and none 
of the developed models have yet won wide acceptance. 
Theoretical attempts to address the splitting failure modes 
of moment-resisting joints do not appear to have been 
reported. 

The present article addresses the splitting failure of the 
simplest possible moment-resisting dowel-type fastener 
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joint, in which the failure is attributed to the perpendicular- 
to-grain loading of one single dowel located close to the 
end of a beam. Beam on elastic foundation (BEF) theory 
for a Timoshenko beam is used, and fracture mechanics is 
introduced through the constitutive relations for the foun¬ 
dation. The stiffness of the foundation is chosen so that the 
perpendicular-to-grain tensile strength and fracture energy 
properties are correctly represented. This particular choice 
of foundation stiffness makes a conventional maximum 
stress failure criterion lead to the same solution as the frac¬ 
ture mechanics compliance method. 


Theory 

Figure 1 shows schematically the type of joint considered, 
its geometry, and assumed failure mode. The encircled part 
in Fig. 1 is assumed to behave as a beam on elastic founda¬ 
tion (BEF), the part above the crack being the beam (edge 
distance, h e , being beam depth), and the part beneath the 
crack serving as the foundation. BEF theory may be ex¬ 
pected to be suitable for application to the problem consid¬ 
ered here because h e in moment-resisting joints is usually 
small as compared with the total beam depth, h. The end 
distance (distance from the critical dowel where splitting 
failure initiates to the beam end) is denoted s. 

The geometry, loading, and boundary conditions for the 
beam on elastic foundation are shown in Fig. 2. Note that P 
is the load on the critical dowel. 

The deflection of the beam axis at the loading point is 
denoted as iv(s'). The foundation stress at the loading point 
is denoted as cr^), and is given by 

ct(s) = Kw(s) (1) 

where K is the foundation stiffness (units: N/m 3 ). 

The deflection at the loading point may in general be 
expressed as 


( 2 ) 








584 



Fig. 1. Definition of geometrical parameters ^*8* 3* Constitutive relation for the foundation 



Fig. 2. Beam on elastic foundation model. Geometry, loading, and 
boundary conditions 
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E being the modulus of elasticity in the grain direction 
and I the moment of inertia of the beam on elastic 
foundation. 

Assuming that the maximum foundation stress occurs at 
the loading point, and that failure occurs when the maxi¬ 
mum foundation stress reaches the perpendicular-to-grain 
tensile strength of the wood, / t , the failure load, P c , thus 
becomes 

FT 1 

P c =f t -( 3 ) 

K W(s) 

The parameters A and 7] are defined as 



Kb 

GA 



Equation 5 is based on the assumption that L —> ©o. This 
assumption does, however, not impose any practical limita¬ 
tions on the model considered here. Previous analyses of 
beams on elastic foundations show that no significant error 
is introduced by assuming L —» oo if E > 10 h e . 


Foundation properties 

The damage and fracture performance of wood is in 
general non linear, but is in the present analysis represented 
by a linear response that is equivalent in terms of peak 
stress, / t , and fracture energy dissipation, ^ f , as indicated in 
Fig. 3. 

The fracture energy of the foundation is in general given 
as 


where b is beam width, I - bh e V^2, A s = 5bhJ6 (rectangular 
cross section, beam depth /z e ), and G is shear modulus. 

The solution to the governing differential equations for a 
Timoshenko beam on elastic foundation divides into two 
solutions depending on A and 77 . By means of solutions 
given by Pilkey, 10 the function *P(s) may for L —» 00 be 
written 


00 

= J odw (6) 

0 

and for the linear response shown in Fig. 3 it follows that the 
stiffness of the foundation, K , is given by 



v 2 =iyiAn, 



The present strength analysis of mode I failure using 
BEF theory and introducing fracture mechanics through 
the linearized constitutive relation as given by Eq. 7 is a 
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complete analogy to the fracture mechanics application of 
the Volkersen model to strength analysis of mode II failure 
in lap joints. 11 The analysis has been termed quasi-non- 
linear fracture mechanics because the material responses 
are assumed to be linear as in linear elastic fracture mecha¬ 
nics (LEFM), but at the same time the tensile strength is 
assigned a finite value, not an infinite value as in LEFM, and 
the finite, nonzero size of the fracture region is considered, 
leading to a failure load, which is not proportional to the 
square root of the fracture energy. It has previously been 
shown 1214 that assigning the foundation the stiffness as 
given by Eq. 7 makes the above-derived solution, which is 
based on a usual maximum stress failure criterion, lead to 
the same solution as obtained by using the compliance 
method of fracture mechanics. 

Two important special cases of Eqs. 3-5 may be identi¬ 
fied, namely s = 0 and 5 —»©©. In both cases, the two solutions 
(Eqs. 5a and 5b) merge into just one solution, which by use 
of Eq. 7 becomes 


Pc = yP< 




2bC l ^[h^, 
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Equation 8 is in agreement with the previously developed 
solutions for glued-in rod joints subjected to pure shear 
load 12 and dowels loaded perpendicular to grain. 14 From Eq. 
8 it follows 



1 c+i P 

2 2£ + l c ’“ 



where P c0 is the failure load for 5 = 0 and P coo is the failure 


load for 5 —»©©. The fact that P c0 < 



for finite values of 


f t is due to the different boundary conditions for the two 
situations; for 5 = 0 the beam end can rotate freely, while for 
s —> °° the rotation of the beam end is zero. 

Results from LEFM are obtained by letting f t —> ©©. In 
Eq. 8, this leads to 


P c , o =^P c ,„=b(GTA/O0 

The LEFM solution obtained here for 5 —> °° is the same as 
that obtained by different LEFM approaches 6,8,15 if the edge 
distance in the latter is small as compared with the total 
beam depth, and it is also the same solution as given by 
another simple LEFM approach 7 if in the latter the shear 
area A s = 5bhJ6 is used instead of the area A = bh e . 



Fig. 4. Influence of perpendicular-to-grain tensile strength (/ t ) on 
failure load (P c ) 
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Fig. 5. Influence of fracture energy (<7 f ) on failure load (P c ) 

In Figs. 4 and 5 are shown examples of the influence of 
perpendicular-to-grain tensile strength, / t , and fracture 
energy, (g t . The examples shown use the following proper¬ 
ties: b = 25 mm, h e - 40 mm, 5 = 160 mm, E = 7200 MPa, G = 
400 MPa; in Fig. 4 Q t = 0.17 N/mm, and in Fig. 5 f t = 1.5 MPa. 
While the failure load asymptotically approaches an upper 
limit for f t —» ©o 5 there is no upper limit for (g t —» ©o. 

It may also be noticed that the present theory predicts a 
very modest influence of the modulus of elasticity (MOE), 
while the shear modulus has a significant impact on the 
failure load. As easily seen from Eq. 8, the failure load 
asymptotically approaches an upper limit for E —> while 
there is no upper limit for G —> ©©. 



Materials and methods 

All specimens were made of glulam of Japanese cedar 
(Cryptomeria japonica ), width, b - 25mm. The laminae 
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Fig. 6. Test setup and geometry of specimens 

(thickness 26 mm) were all the same grade and were without 
finger joints. MOE in the grain direction was determined 
by measuring the longitudinal vibration frequency of the 
glulam beams, from which the specimens were cut. All tests 
were displacement controlled and failure occurred after 
2-3 min. 

The symmetrical three-member joints, with the glulam 
test specimen as the middle member and 25-mm plywood 
plates as side members, were subjected to a pure moment in 
a four-point bending test setup as shown in Fig. 6, causing 
splitting failure at the rightmost dowel (edge distance h e , 
and end distance s). Dowels with diameters of 14 mm were 
used in 15-mm holes, and a gap was deliberately made 
between middle and side members to assure no transfer of 
forces due to friction. 

Nine different combinations of h e and s were tested: For 
h e = 20 mm and for h e = 40 mm, s = 20, 40, 80, and 160 mm 
were tested, and for h e - 60 mm, s = 40 mm was tested. 
Twelve specimens were prepared for each condition. How¬ 
ever, only specimens without knots in the vicinity of the 
point of crack initiation were selected in order to limit the 
variation, resulting in 7-12 specimens tested for each condi¬ 
tion. Moisture content (MC) was 14%, density was 
372 kg/m 1 at the given MC, and MOE was 7200 MPa. 

The load, P, acting perpendicular to grain on the critical 
dowel is calculated as 

p = 2AP l + p : ; w (10) 

where 2 P L is the total load applied by the testing machine 
and P sw ~ 160 N is the contribution from the weight of the 
plywood. 

Tests were also conducted on double cantilever beam 
(DCB) specimens (Fig. 7a) to determine the fracture energy 
of the glulam, and on plate joint specimens (Fig. 7b) to 
determine the tensile strength. Eleven knot-free DCB 
specimens, b = 25 mm, h - 40 mm, a = 250 mm, and L - 
550mm, were tested. Plate joints were tested using two 
different edge distances: h e = 20 mm and 40 mm. Ten to 
12 knot-free specimens with b - 25 mm, h - 210 mm, and 
L - 500 mm were tested for each h e value. 


Results and discussion 

The model derived in the present article requires two 
fracture properties, fracture energy (g f and perpendicular- 



Fig. 7. Double cantilever beam specimen for determination of fracture 
energy (a) and plate joint specimen for determination of perpendicu- 
lar-to-grain tensile strength (b) 


to-grain tensile strength / t , as input parameters. Because f t 
is highly volume dependent, it is not obvious what volume 
(or cross-sectional area) should be used if attempting 
to determine it directly by tension tests. In the following, f t 
is therefore determined by means of the plate joint 
specimen. 

Analyses of the DCB and plate joint specimens based on 
the same BEF and fracture mechanics theory as given in 
this article have been presented elsewhere. 13,14 These mod¬ 
els likewise both depend on (g f and / t , and an iterative pro¬ 
cess of successive calculations is necessary (two to three 
calculations) to derive (g f and f t . The results obtained were 
( g f = 0.21 N/mm and f t = 1.05 MPa, where f t is the mean value 
obtained from plate joint specimens with h e = 20mm (/ t = 
1.29MPa) and 40mm (/ t = 0.81 MPa). 

The fracture energy may also be determined from the 
DCB specimen tests using a simple linear elastic fracture 
mechanics solution, which does not depend on / t , resulting 
in (g t = 0.17 N/mm. In that case the plate joint tests 14 result in 
f t = 1.64 MPa (mean of 1.22 MPa and 2.05 MPa obtained for 
h e = 40mm and h e - 20mm, respectively). In the following 
(/ t , (g f ) = (1.05 MPa, 0.21 N/mm) is used, but (/ t , (g f ) = 
(1.64 MPa, 0.17 N/mm) also leads to good predictions. The 
fact that significantly different values of f t are obtained for 
h e = 20 mm and 40 mm questions the suitability of using the 
plate joint specimen for deriving / t . Further discussions on 
this problem may be found elsewhere. 13,14 

Figure 8 shows the results of the tests on moment- 
resisting joints compared with the predictions as given by 
Eq. 5a. The theoretical model predicts amazingly well the 
influence of edge distance, h e , and end distance, s. 

Figure 8 indicates that a bilinear approximation may be 
suitable for a simplification of Eq. 5a. Equation 8 gives the 
point for s = 0 and the horizontal upper limit. Using the 
tangent at s = 0 leads to the simple expression 
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Fig. 8. Equation 5a compared with experimental failure loads. f = 
1.05 MPa, (g t = 0.21 N/mm, bars indicate standard deviation 



Fig. 9. Equation 11a compared with experimental failure loads. 
f t = 1.05 MPa, (g { = 0.21 N/mm, bars indicate standard deviation 
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Equation 11a is shown together with the test results in Fig. 
9. Judging from Fig. 9, the bilinear relation given by Eq. 11a 
seems in general to give a slightly conservative (as com¬ 
pared with Eq. 5a) but fairly good approximation, and may 
be worth considering as a simple design tool. While the 


solutions given by Eqs. 5a and 5b are each numerically 
restricted to cover certain combinations of some geometri¬ 
cal and material properties, Eq. 11a contains no such re¬ 
strictions and may, without numerical problems, be used for 
any finite geometrical and material properties. 

Although Eq. 11a leads to very good predictions for 
the joints reported in the present article, it still remains 
to be investigated whether it holds for other conditions as 
well. 


Conclusions 

A quasi-non-linear fracture mechanics model based on 
BEF theory was applied for the analysis of splitting failure 
in moment-resisting dowel joints. An approximation simple 
enough for practical design was also proposed. 

Tests were conducted on simple moment-resisting 
joints with only two dowels, and end distance and edge 
distance of the critical dowel were varied. Although theo¬ 
retical predictions were in excellent agreement with test 
results, indicating that the theoretical approach is basically 
sound, further testing is needed on more realistic dowel 
configurations. 

The models presented depend on the perpendicular-to- 
grain tensile strength. The tensile strength is highly volume 
dependent, and it is difficult to determine it directly by 
tension tests because the appropriate volume (or cross- 
sectional area) of the test specimen is not obvious. In the 
present article, so-called plate joint specimens were used 
for estimation of the tensile strength. More investigations 
are needed to reveal whether the plate joint specimen 
is suitable as a standard test specimen for that purpose 
or if another way of deriving an appropriate value of the 
perpendicular-to-grain tensile strength is needed. 
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